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ABSTRACT 

We report the discovery of a pair of quasars at z = 1.487, with a separation of 8'.' 585 ± 0'.'002. 
Subaru Telescope infrared imaging reveals the presence of an elliptical and a disk-like galaxy located 
almost symmetrically between the quasars, in a cross-like configuration. Based on absorption lines 
in the quasar spectra and the colors of the galaxies, we estimate that both galaxies are located at 
rcdshift z = 0.899. This, as well as the similarity of the quasar spectra, suggests that the system is 
a single quasar multiply imaged by a galaxy group or cluster acting as a gravitational lens, although 
the possibility of a binary quasar cannot be fully excluded. We show that the gravitational lensing 
hypothesis implies that these galaxies are not isolated, but must be embedded in a dark matter halo of 
virial mass ~ 4 x 10 14 Vq 1 Mq assuming an NFW model with a concentration parameter of c V i r = 6, or 
a singular isothermal sphere profile with a velocity dispersion of ~ 670 km s^ 1 . We place constraints 
on the location of the dark matter halo, as well as the velocity dispersions of the galaxies. In addition, 
we discuss the influence of differential reddening, microlensing and intrinsic variability on the quasar 
spectra and broadband photometry. 

Subject headings: gravitational lensing: strong — quasars: individual: SDSS J132059. 17+164402. 59 
- quasars: individual: SDSS J132059.73+164405.6 



1. INTRODUCTION 

Quasar pairs belong to one of three categories: gravi- 
tational lenses, binary quasars, and apparent (projected) 
pairs. Both components of a pair in the first two 
categories have the same redshiftP^l. Gravitationally 
lensed quasars dominate the small separation A9 < 3" 
quasar pairs. These are produced by galaxy scale lenses, 
and for the past three decades have proven to be in- 
valuable astrophysical and cosmological probes (e.g., 
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of the UH2.2 m telescope for the observations is supported by 
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12 Small redshift differences in the components of binary quasar 
pairs such as 661 ±173 km s" 1 for LBQS 0015+0239 (Impcy ct al. 
120021 ). are consistent with the line-of-sight velocity difference of 



IClaeskens fc Surdeil I2002t ISchneider etaLl [2006). At 
larger separations, the lensing probability falls quickly 
((Turner et al.l I1984D . and becomes increasingly domi- 
nated by envi ronmental effects du e to galaxy groups and 
clusters (e.g.. iKeeton et all 120001) . The largest separa- 
tion, double-image lens ed quasar confirm ed to date is 
Q0957+561 (A0 = 6"17. IWalsh eFa!lT979l) . and the only 
larger separation lensed quasars known are the five - image 
SDSS J1004+4112 (A0 = 14'/ 7. Ilnada etld1l2Tj(jl l2005h 
and the three-image SDSS J1029+2 623 (A0 = 22'.'5, 
Ilnada et al.l 120061: IQguri et al.ll2008bl ). The two latter 
systems are lensed by clusters. Clearly, a larger sam- 
ple is needed to study observationally the lensing prob- 
ability distribution in the large-separation range (e.g., 
iQguril 120061: iMore et al.l 120121) . These wide lenses pro- 
vide unique opportunities to study the interplay between 
baryonic and dark matter in groups (AO < 10"), as 
well as the dark matter distribution in clusterj^l (e.g., 
IQguri et alJl200l . 

Most large-separation quasar pairs known, however, 
are binary quas ars (e.g., 19 pairs wit h 3" < A(9 < 10" 
in the sample of lHennawi et al.l 120061 ). These are physi- 
cally associated quasars, in the sense that they are either 
gravitationally interacting, or belong to the same group 
or cluster. Binary quasars are found to be more abun- 
dant than predicted from extrapolatin g the quasar cor- 
relation function to small scales (e.g.. iDiorgovskil Il99ll: 
Hewett et all 119981 : iHennawi et "all 120061 : iKavo fc Oguril 



bound pairs of galaxies. 



20121) . which may imply that quasar activity is trig- 
gered and susta i ned by tidal interactions or merger s (e.g. , 
I Mortlock et all fl999t IHopkins et al.1 [20051 1200a. [20071 
2008: IGreen et al.H2010L and references therein). 
Many quasar pairs found in the range 3" < A9 < 10" 

13 It is worth mentioning that, in addition to lensed quasars, 
there are also many known lensed galaxies, with a lar ge range of 
image separations, which we do not addre ss here fe.g,. lFaure et all 
120081: lAuger et afll200Bl: IMore et al.ll2012Tl . 
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have proven difficult to classify as either lensed or b i- 
nary quasars (e.g. Q2345+007; IWeedman et~aTl [1982) . 
In these systems no plausible lenses have been identified 
in the foreground, although the redshifts of the two com- 
ponents are identical and the spectra are quite similar. 
The spectra of gravitationally lensed images need not 
necessarily be identical, since they are prone to effects 
such as delayed intrinsic variability, chromatic microlens- 
ing, and differen tial extinction in the le nsing g alaxy (e.g., 
iWucknitz et~all 120031 : lYonehara et~aT1 120081 ). On the 
other hand, there is little dynamic range in the spectral 
variation of quasars, me aning that pairs can h ave similar 
spectra by chance (e.g., Mortloc k et al.|[l999h . 

The suggestio n of lensing by "dark" lenses (e.g., 
iKoopman s 2000) or cosmic strings (|VilenkinlH98l has 
been put forward for these quasar pairs. There are in 
fact very few criteria for discriminating between the lens- 
ing and binary quasar hypotheses, without monitoring 
the luminosity variability of the components, in order to 
look for a time delay. Instead, clues about this popula- 
tion of pairs have been found statistically. For instance, 
it has been argued that most of these quasar pairs are 
not lensed quasars, based on compariso ns of the optical 
and r adio properties of the population (Kochanc k et al.l 
1999). As for individual pairs, their nature must in gen- 
eral be estimated from the ensemble of their observed 
properties. 

In this paper, we report the discovery in the 
course of the Sloa n Digital Sky S urvey Quasar 
Lens Search /SOLS : IQeuri et all [200l [2008al I2012bt 
llnada et all [20081 120101 I2012T ) of a new large- 
separation (A8 ~ 8'.'5) quasar pair which has 
proven difficult to classify, SDSS J132059. 17+164402. 59 
and SDSS J132059. 73+164405.6 (hereafter jointly 
SDSS J1320+1644). In the near infrared, we detect two 
galaxies located between the quasars, producing an ap- 
parent cross-like configuration. In <J21 we give a brief de- 
scription of the lens candidate selection from the SDSS 
data, and in <j3]we present our imaging and spectroscopic 
follow-up observations. Next, we estimate photometric 
redshifts for the two galaxies and show that they are lo- 
cated in the foreground (Sg]). We subsequently proceed 
with the gravitational lens mass modeling of this system 
(§SJ). In fJHl we examine the quasar spectral and photo- 
metric differences, whereas in ^7|we explore the environ- 
ment of the system. Finally, we summarize our conclu- 
sions in SjSJ Throughout this paper, we assume the con- 
cordance cosmology with Hq = 70 h-?o km -1 s _1 Mpc -1 , 
M = 0.27 and n A = 0.73. 

2. DISCOVERY IN THE SLOAN DIGITAL SKY SURVEY 

The Sloan Digit al Sky Survey ( SDSS-I, 2000-2005, 
SDSS-II, 2005-2008;|York_etal. 2000) is a combination of 
imaging and spectroscopic surveys which have mapped 
~ 10,000 square degrees of the sky, centered at the 
North Galactic Cap. The observations are made at the 
Apache Point Observatory in New Mexico, USA, with 
a ded icated 2.5-meter wide-field telescope (|Gunn et all 
2006). The imaging survey is conducted in five broad- 
band filters, centered at 3561(m), 4676 (g), 6176 (r), 7494 
(z), a nd 8873A (z) (i Fukugita et all 119961: iGunn et all 
119981: iDoi et alll2010l) . The automated pipeline reduc- 
tion achieves an astro metric accuracy better than about 
O'.'l (|Pier et al.l 120031 ). and a photometric zeropoint ac- 



curacy better than about 0.01 magnitude o ver the en- 
tire survey area, in the g, r, and i bands dHogg et al 
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2006; 



ISmith et al.ll2002Ulvezic et al.l l2004al iTicker et al 
IPadm anabh an et al.l I2008D . The spectroscopic 
observations are carried out with a multi-fiber spec- 
trograph between 38 00 — 9200 A, with a resolution of 
R ~ 1800 - 2100 (|Blanton et all l200l . All data 
have bee n made publicly avai l able, in periodic data 
relea s es (IStoughton et al.l 120021: lAbazaiian et all 1 2003 . 
2001 120051 120091: lAdelman-McCarthv et alll2006l 12007. 
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The SQLS has been one of the most successful strong 
lens surveys conducted to date. It identifies lensed quasar 
candidates among spectroscopically confirmed quasars i n 
the SDSS (|Richards et al.ll200l ISchneider et alJl20loft . 
by combining two selection algorithms, designed to find 
both small- and large-separation lenses, respectively. 
The morphological selection identifies small-separation 
candidates that are unresolved, yet poorly fitted by the 
point spread function (PSF); the color selection algo- 
rithm searches, in the vicinity of every confirmed quasar, 
for resolved stellar objects that have similar colors to 
the quasar. The SQLS has discovered m ore than 40 new 
lensed quasars so far (jlnada et al.ll2012T ). Together with 
previously known lenses rediscovered by the SQLS, it has 
produced a sample of ~ 60 lensed quasars, roughly half 
of all lensed quasars discovered to date. 

SDSS J1320+1644 is a large-separation lensed quasar 
candidate identified by the color selection algorithm. We 
reproduce the SDSS multi-band composite image of the 
system in Figure [TJ A faint red object, classified by the 
automatic pipeline as a galaxy, is detected in between 
the two blue stellar components constituting the lensed 
quasar candidate, towards North. We give the SDSS pho- 
tometric results for the quasars and the red object in Ta- 
ble [1] As a result of the finite size of the optical fibers in 
the SDSS spectrograph, spectra are usually not obtained 
for two objects less than 55" apart (the "fiber collision"; 
IBlanton et al.ll2003l ). and therefore SDSS spectroscopy is 
available for only one member of the pair. We repro- 
duce the spectrum of that object, classified as a quasar 
at z = 1.5024 ± 0.0024, in Figure [2] (in blue). 

3. FOLLOW-UP SPECTROSCOPY AND IMAGING 

3.1. Spectroscopy of the quasar pair 

Follow-up spectroscopy aimed to ensure that both stel- 
lar components are quasars at the same redshift, and to 
assess their spectral similarities, was performed at the 
Astrophysical Research Consortium Telescope (ARC 3.5 
m), located at the Apache Point Observatory. Spec- 
tra were taken with the Dual Imaging Spectrograph, 
equipped with a combination of blue and red gratings, 
B400/R300. The observations were conducted on 2009 
February 19 with a 1'.'5 width slit oriented so that both 
stellar components are on the slit. One 1500 s exposure 
was taken at an effective airmass of 1.63. The extracted 
spectral coverage was 3700 - 10000 A, with a resolution 
R ~ 500. The data reduction was performed using stan- 
dard IRAF0 tasks. 

14 The interactive Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, which 
are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National 
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The spectra shown in Figure [5] (black and red) indi- 
cate that both stellar components have similarly shaped 
quasar broad emission lines (BELs; C III], Mg II) at the 
same wavelengths, and therefore are quasars at the same 
rcdshift. Our redshift estimate based on the C III] and 
Mg II peaks (including a reanalysis of the original SDSS 
spectrum, which was taken in a longer, 4200s exposure, 
and better resolution R - 1850 - 2000) is 1.487±0.001, 
which is slightly lower than, and incompatible with the 
value in the SDSS database and its associated error bar. 
We identify absorption lines in both quasars as the Mg II 
2796 A, 2803 A doublet and possibly Fe II 2599.4 A, in- 
dicating z = 0.899. We identify another Mg II doublet 
at z = 1.168, present only in quasar B (the quasar des- 
ignation corresponds to that in Figure 2]) . 

The flux ratio A/B increases with wavelength, from ~ 
0.8 to ~ 1.5. This could be explained by a combina- 
tion of microlensing, differential extinction and intrinsic 
variability, or simply differential slit losses, and does not 
constitute by itself an argument against the gravitational 
lensing hypothesis. We provide a more thorough compar- 
ative analysis of the spectra in S£SJ 

Overall, the identical redshift of the two quasars, the 
similarity of the shapes of the BELs, as well as the pres- 
ence of absorption lines at the same redshift, appear 
consistent with the gravitational lensing hypothesis. Of 
course, these features are fully consistent with the binary 
quasar hypothesis as well, in which case spectral differ- 
ences arise naturally, as there are two distinct quasars. 

3.2. UH88 and Subaru Telescope imaging data 

Optical follow-up imaging observations were originally 
conducted with the Tektronix 2048x2048 CCD camera 
(Tek2k; V, R, I and z bands) and the Wide Field Grism 
Spectrograph 2 (WFGS2; z band) at the University of 
Hawaii 2.2-meter (UH88) telescope, in April 2009. The 
Tek2k camera has a 7'. 5 x 7'. 5 field of view and a pixel 
scale of 07219 pixel" 1 . The WFGS2 has a 11' x 11' field 
of view and focal reducer of 0734 pixel -1 , and was used 
in imaging mode with the Tek2k camera. Bias and flat 
frames were also obtained, and the data were reduced 
using standar d IRAF tasks . The standard star system 
PG0918+029 (|Landoltlll99l ) was observed for photomet- 
ric calibration, but the zero-points determined from the 
five individual stars have a scatter of ~ 0.2 mag, possi- 
bly indicative of non-photometric conditions. In order to 
obtain a better precision, the instrumental magnitudes 
of 10 bright stars in the field of view were compared to 
their SDSS ugriz magnitudes. The Lupton (2005) for- 
mulaJ^I were used to transform between ugriz and VRI 
magnitudes, which resulted in smaller scatter. The ex- 
posure time, airmass, seeing, zero-point uncertainty and 
observation date are given for each band in Table [2j 

In the deeper UH88 images (the I and z filters), in 
addition to the faint photometric object detected be- 
tween the two quasars in the SDSS data, another ob- 
ject is detected, located almost symmetrically with re- 
spect to the quasars (Figure [3]). In order to determine 
the morphological and photometric properties of these 
two objects, as well as to estimate photometric redshifts, 
we obtained higher resolution near-infrared (J, H and 

Science Foundation. 

15 http : //www . sdss . org/dr5/algorithms/sdssUBVRITransf orm .html 



Ks bands) images with the Multi- Object InfraRed Cam- 
era and Spectrogr aph (MOIRCS; llchikawa et~all 120061: 
iSuzuki et al.l [20 08) . at the Subaru 8.2-meter telescope 
(jive et a 1.1 1200 Hi. in April 2010. MOIRCS has two 
HAWAII-2 2048x2048 detectors providing a 3'.94 x 6'. 90 
total field of view, with a pixel scale of 07117 pixel -1 . 
The data (Tabled]) were reduced with the MCSRED soft- 
ware package (I. Ta naka et al., in prepa ration), and the 
standard star FS33 ([Leggett et al.ll2006h was used to es- 
timate the photometric zero points. 

In the J, H and Ks imaging (Figure [4}, the two faint 
objects are clearly detected. They have extended mor- 
phology and red colors, indicative of fairly high-redshift 
galaxies. In each band, we modeled all four ob jects simul- 
taneo usly using the public software GALFIT ([Peng et al.l 
2002), with nearby bright stars as PSF templates. We 
fitted the galaxies with a Sersic profile convolved with the 
PSF. The resulting astrometry and photometry (includ- 
ing the UH88 photometry) are given in Tables [3] and [H 
respectively. We were able to characterize the morphol- 
ogy of the two galaxies, which we summarize in Table 
[H The Sersic index of Gl is close to the canonical value 
for elliptical galaxies, n = 4, whereas for G2 it is close 
to the typical value for disk galaxies, n = 1. We there- 
fore assume that Gl is an elliptical galaxy, whereas G2 
is disk-like. This is also consistent with the bluer color 
of G2. 

We also detected three additional fainter objects in the 
vicinity of the system. While FWHM measurements and 
GALFIT modeling are not very reliable, our analysis sug- 
gests that the objects are extended, and so we mark them 
as galaxies G3, G4 and G5 in Figure [U Aperture pho- 
tometry, summarized in Tabled indicates that the colors 
of G4 (but not of G3 and G5) are very similar to those 
of Gl and G2, and we therefore assume that this is a 
galaxy at about the same redshift. We also remark that 
the J — H colors of G3 and G5 are different from those 
of A and B, and it is therefore unlikely that these are 
additional, faint gravitationally lensed quasar images. 

In the UH88 and Subaru Telescope imaging observa- 
tions, the flux ratio A/B is above unity, in opposition to 
the SDSS results. The chromatic change in the bright- 
ness ranking of the quasar images could be explained by 
microlensing (see 



4. PHOTOMETRIC REDSHIFTS 

In order to determine whether the two brighter galax- 
ies are located in the foreground of the quasar pair, 
as is required by the gravitational lens hypothesis, we 
need to estimate their redshifts. We estimate pho- 
tometric redshifts, based on the magnitudes of the 
galaxies in different filters, using the template-fitting 
meth ods implemented in the publicly available Hy- 
perZ ([Bolzonella et al.l l2000) and EAzY ([Brammer et al.l 
2008) algorithms. These algorithms fit the observed mag- 
nitudes to spectral energy distribution (SED) templates, 
via x 2 minimization. As templates, we employ the ob- 
served mean SEDs of local galaxies from iColeman et al.l 
(1980, hereafter CWW), which are extrapolated into the 
ultrav iolet and near-infrared wit h the evolutionary mod- 
els of lBruzual fc Chariot! (| 19931 ) . The EAzY algorithm 
also uses the redshift distribution of galaxies of a given 
appare nt magn i tude a s a Bayesian luminosity prior, fol- 
lowing iBenitezI (|2000l) . This helps to break the degen- 
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eracies between multiple probability distribution peaks 
at different redshifts. 

In order to obtain consistent results for the two al- 
gorithms, we had to increase the uncertainties on the 
magnitudes to about double the values quoted in Table 
[4j which also makes most of the aperture- and model- 
determined values consistent. In addition to the statis- 
tical errors quoted in that table, additional sources of 
errors could potentially be introduced by the zero-point 
uncertainties (Table [2J , possible blending of sources due 
to their close proximity and large seeing, different pixel 
size in the V, R, I and z band observations (aperture 
photometry), or the modeling difficulties and PSF un- 
certainties (model magnitudes). 

The results of the photometric redshift estimates are 
shown in Figure [5] and Table El Gl is well- fitted by an el- 
liptical template, and G2 by an Sbc template, consistent 
with the expectations from the galaxy morphologies (see 
£|3.2j) . Although the probability distributions are rather 
broad, they are in agreement, inside the la interval es- 
timated from the redshift probability distributions, with 
z = 0.899, which is the redshift of the absorption lines 
identified in both spectra. The redshift of the second 
Mg II absorption system in quasar B, at z = 1.166, is 
rejected at about 2a confidence level, while the proba- 
bility of a redshift as high as the quasar redshift appears 
negligible. Therefore, throughout the rest of this work, 
we assume that these galaxies are responsible for the ab- 
sorption lines at z = 0.899. 

For the faintest objects we detected, G3, G4 and G5, 
the redshift probability distributions are very broad (Ta- 
ble [5]) . G3 and G5 show a smaller preferred redshift of 
z ~ 0.5. On the other hand, G4 has a preferred redshift 
of 1 . 1 1 , very close to that of the Mg II absorption system 
at z = 1.168. This absorption system is quite strong, 
with equivalent width ~ 2.9 A in the rest frame of the 
absorber. Here the projected proper distance between 
B and G4 is 21 h? Q kpc at z = 1.168. On the other 
hand, the Mg II absorption system at z = 0.899 has 
equivalent widths ~ 2.9 A and ~ 4.5 A (in the spectra 
of quasars A and B, respectively) for impact parameter 
~ 35 /lyp 1 kpc, corresponding to the distances between 
the quasars and Gl, G2. Such strong absorptions at 
small imp act parameters are co nsistent with the litera- 
ture (e.g., iChurchill et al.ll2005l and references therein). 
The stronger absorption in quasar B could be explained 
by its closer proximity to one of the galaxies (Gl). 

5. GRAVITATIONAL LENS MASS MODELING 

We next proceed with modeling the system as a gravi- 
tational lens. The observational constraints provided by 
the available astrometry and photometry are the posi- 
tions of the images A and B, their fluxes (we take the 
flux ratio to be either ~ 1.0 ± 0.2 or ~ 1.4 ± 0.2, as 
we will show in and the positions of the two main 
galaxies expected to act as lenses, Gl and G2. We also 
know the redshift of the source, z — 1.487, and of the 
lenses, z = 0.899. As the number of observational con- 
straints is limited, we focus on simple SIS and NFW lens 
models and their elliptical counterparts, with or without 
shear. The singular isothermal profile is known to be a 
good approximation for the mas s distribution in galaxies 
(e.g., iRusin fc Kochanekl 120051 and references therein). 



It has also been employed in the le ns modeling of galaxy 
clusters (e.g., Ilnada et all l200l. The NFW profile 
(N avarro on the other hand, is frequently 

employed as an approximation of the mass distribution 
in dark matter haloes, which dominate galaxy clusters. 
Constraints can also be placed on the "strength" of the 
lenses, represented (in the case of the SIS model) by their 
velocity dispersions, which we infer from the observed 
luminosities. To this end, we estimate the rest frame R 
band absolute magnitudes of Gl and G2 as ~ —22.5 and 
~ —22.3, respectively, calculated with HyperZ from the 
best- fit E and Sbc templates . We use the Faber- Jackson 
(F -J; iFaber fc Jackson! [19761) and the Tully-Fisher (T- 



F; iTullv fc FisheJ 19771) la ws, with parameters given in 
iBinnev fc Tremaind ( 20081 ) , and infer velocity dispersions 
of 237 ± 54 km s" 1 and 163 ± 15 km s _1 , respectively. 
For G2, which is disk-like, we have assumed a singular 
isothermal sphere (SIS) profile and estimated the veloc- 
ity dispersion as a = Wcirc/V^j where w c ; rc is the circular 
velocity obtained from the T-F law. 

The velocity dispersion estimate above implies that the 
brightest galaxy Gl has an Einstein radius of only ~ 
/ .'5±0 / /2, and an inclosed mass of - (2.1±1.5) x 10 11 
Mq. However, the large separation between A and B 
(~ 8'.'6), assuming they are lensed images of a single 
source, implies an enclosed mass of - 1.2 xlO 13 M . 
Consequently, this cannot be a gravitational lens system 
where the main deflectors are the individual members of 
a group of galaxies, and we need to invoke a significant 
embedding dark matter halo. 

A question worth posing is to what extent must the 
faintest objects G3, G4 and G5, be taken into account 
when building gravitational lens models. We only es- 
timate the velocity dispersion of G4, assumed to be at 
z = 0.899 orz= 1.168 (see gp. We neglect the contri- 
bution of G3 and G5, as they are located further away 
from images A and B, and are most likely at a signifi- 
cantly different redshift from Gl, G2 and G4. Assuming 
that G4 is located at the same redshift as Gl and G2 
leads to velocity dispersions a — 150 ± 35 km s _1 or 
a = 105 ± 10 km s _1 , from the F-J and T-F laws, re- 
spectively. A smaller effective velocity dispersion w ould 
be inferred if placing G4 at redshift z = 1.168 (see E1A.1I 
for details). 

We subsequently explore gravitational lens models in 
which the lensing potentials of Gl, G2 and G4 are 
boosted by an embedding dark matter halo. We model 
the observed galaxies with SIS profiles, and the dark mat- 
ter halo with either an SIS or an NFW profile. We con- 
sider the case where the dark matter halo is centered at 
the location of the brightest galaxy Gl, as well as the 
case where its position is a free parameter, on which we 
obtain constraints. We only consider models with de- 
grees of freedom (d.o.f.), which we are able to fit perfectly 
(x 2 € 1). This means that for the models where the cen- 
ter of the dark matter halo is fixed, we save two parame- 
ters which we can use to introduce a shear / ellipticity and 
associated position angle. T he gravitatio nal lens model- 
ing is performed with glafic ( Oguri 2010). 

We show the critical lines and caustics for these mod- 
els in Figure [6] We report the best fitted parameters of 
the models (velocity dispersion or virial mass of the dark 
matter halo, ellipticity/shear and associated position an- 
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gle, total magnification and time delay) in Table [7J The 
velocity dispersions of G2 and G4 (also Gl in the case of 
the "SIS free" and "NFW free" models of Figure EJ were 
fixed at the values inferred earlier, but we also explored, 
on a suitable grid, combinations of velocity dispersions in 
the ranges allowed by the F-J and T-F laws. In this case, 
we only considered models with x 2 < 2.3, corresponding 
to the la confidence interval for two parameters. The 
effect of the range of velocity dispersions is reflected in 
the range of values for the parameters in Table [71 as well 
as the allowed position of the dark matter halo, for the 
'"SIS free" and "NFW free" models. 

There are two major differences between the NFW and 
SIS models considered for the dark matter halo. First, we 
parameterize the NFW profile by its virial mass M„j r and 
concentration parameter c„j r (as opposed to just a sin- 
gle parameter, a, for the SIS model), and we fix c V i r at a 
fiducial value 6, close to the typical value for galaxy clus- 
ters a t the mass range and reds hift we are interested in 
(e.g., iBinnev fc Tr emaine 2008, and r efere nce therein). 
We explore different choices of c V i r in Second, the 

NFW profile is shallower than SIS at the center, and 
can therefore produce an additional faint central image. 
Therefore, for the "NFW free" model, we require that 
the model produce no more than three images, and that 
the central image be ~ 50 times or more fainter than B, 
corresponding to our detection limit at S/N ~ 5 in the 
deepest band, Ks. The central image cannot be one of 
the faint objects G3 or G5, as those objects are too far 
away from the midpoint between A and B, are also most 
likely extended, and have different colors. 

Finally, in addition to the models where the dark mat- 
ter halo is fixed on Gl, we tried a similar range of SIS+7 
models at the location of G2, but were unable to produce 
good fits. We note that G2 is farther away from the line 
connecting images A and B, and therefore a model cen- 
tered on G2 is expected to do a poorer job of fitting the 
data. We do not try to fix the dark matter halo at the 
location of G4, as we are not sure whether this galaxy is 
at the same redshift as Gl and G2, and we also do not 
expect the dark matter halo to be centered on a fainter 
galaxy of the group/cluster. 

Our main result in this section is that any reasonable 
mass model reproduces the observed image configura- 
tion. In order for this system to be a two-image gravita- 
tionally lensed quasar, we do however require a massive 
dark matter halo in which the visible galaxies are em- 
bedded. This halo should be centered either on Gl (but 
not on G2) or at some location close to the midpoint 
between A and B. The choice of the intrinsic flux ratio 
(either ~ 1.0 or ~ 1.4) does not a ffect t he successful mass 
models qualitatively. Finally, in ^A.3I we discuss several 
caveats regarding the estimated velocity dispersions of 
the galaxies. 

6. QUASAR SPECTRO-PHOTOMETRIC DIFFERENCES 

The three physical phenomena known to modify the 
spectra of the images of gravitationally lensed quasars are 
differ ential extinction due to dust in the lensing galaxy 
(e.g., iFalco et aLlll999D, microlensing due to stars in th e 
lens (e.g.. lChang fc RefedaHU979t ISchneider et al.ll2006l) . 
and intrinsic quasar variability coupled with the time 
delays between images. The first process is a static one, 
whereas the other two are time-dependent. We assess 



their effects on the continuum and broad-band fluxes in 
turn. 

6.1. Differential extinction 

The most striking difference in the follow-up spectra 
of quasars A and B is the gradual drop in the flux ra- 
tio A/B with wavelength, shortward of about 7000 A. 
Such an effect has previously been observe d in several 
lense d quasars, such as SDSS J 1001+5027 (IQguri et al.1 
l2005h and SDSS J1313+5151 (IQfek et all 120071) . Dif- 
ferential extinction (reddening) between the quasar im- 
ages could explain this effect, as it is most effective at 
short wavelengths. Here we caution however that the 
monotonous drop in the flux ratio may also be due to 
unknown differential slit losses in the ARC 3.5 m spec- 
troscopy. Nonetheless, we do not expect the differential 
slit losses to be large, because the flux drop is supported 
by the available photometry from u to z band, which 
shows quasar A to be redder than B (see also Figures 
Q] and [JJ . In addition, the flux ratios inferred from the 
UH88 photometry, obtained just two months after the 
follow-up spectra, are consistent with the spectral ra- 
tios at the respective wavelengths. On the other hand, 
extinction should be weakest in the longest wavelength 
filter Ks, which corresponds roughly to J band in the 
rest frame of Gl and G2. A caveat however is that there 
is no obvious galaxy responsible for the extinction in A, 
i.e. located much closer to A than the other galaxies are 
to B. 

In contrast, the large change in the spectral slope of 
quasar A at short wavelengths, on a time span of 13 
months, between the original SDSS and ARC observa- 
tions, requires a significant contribution of the time- 
dependent phenomena. While the slope change can also 
be caused by slit losses, the two sets of photometric data 
in the V - z band range, taken ~ 4 years apart, also 
shows the brightness rank of the two quasars interchang- 
ing, with A getting brighter, and B becoming fainter. 
Therefore time-dependent phenomena must be acting on 
the quasars. We next look into the extent to which mi- 
crolensing can explain these data. 

6.2. Microlensing 

We assume that the microlensing variability is due to 
the motion of stars in the lensing galaxies, which we 
consider to be dominated by the proper motion of the 
galaxies, with transverse velocity v± ~ a, where a is 
the vel ocity dispersion in the galaxy group/cluster. Fol- 
lowing Richar ds et al.l (|2004f ) in the following formulae, 
there are two timescales to the microlensing event. The 
first is the expected duration of the event, which is the 
time in which the microlens crosses the source, 

where i? src is the size of the source quasar accretion disk 
in the rest frame ultraviolet, corresponding to the optical 
region of the observed spectrum of SDSS J1320+1644. It 
has an estimated ex tent of a few light-days or less (e.g., 
iMorgan et all [20081 ). 

The second timescale is the mean interval between mi- 
crolensing events, which is crudely estimated as the time 
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needed to cross an Einstein radius, 

/ M \ 1/2 / v x \ _1 1/2 

(2) 

with M being the mass of the microlens. 

Although these results are rough estimates, they are 
close to the timescales on which our data shows spectral 
and photometric variations. Careful analytical consid- 
erations of the expec ted microlensing light curves (e.g., 
lYonehara et al.ll2008l ) as well as ray-shooting techniques 
to account for the fact that the source could be mi- 
crole nsed by an ensembl e of stars at any particular time 
(e.g.. iKavser et aLl fl986) show that microlensing can in- 
deed cause variations in amplitude comparable to, or 
even larger than we observe in the spectro-ph otometric 
data. In standard accretion disk theory (e.g., iPetersonl 
Il997f ). short wavelength emission originates at compar- 
atively smaller distances from the central engine, and 
so is preferentially affected by microlensing. This is in 
agreement with the large changes we detect at progres- 
sively shorter wavelengths in our two-epoch spectra of 
quasar A. Conversely, since the near-infrared emission is 
expected to originate at larger distances from the cen- 
ter than the UV or optical continuum, we expect the 
longest wavelength band Ks (originating approximately 
at z band in the quasar rest frame) to be less affected by 
microlensing. 

On the other hand, the BEL regi on is as large as a 
few light-weeks or light- months (e.g.. lClavel et al.lll991h . 
with lower ionization lines originating at larger distances. 
Thus, we expect that, in particular, the Mg II BEL is less 
prone to microlensing effects. We find that the spectral 
flux ratio A/B of the two images (Figure [2]) shows a de- 
crease at the center of the Mg II BEL, compared to the 
surrounding continuurrF^I. This could be indicative of 
microlensing affecting the continuum, but not the Mg II 
BEL. 

In principle, an estimate of the intrinsic flux ratio 
can be obtained from the flux ratio in the continuum- 
subtracted BELs, since they are less prone to microlens- 
ing effects than the continuum fluxes. However, in our 
case there is the chance that differential slit losses, al- 
though believed to be small (see £16.11) affect the BELs, 
and therefore we refrain from doing so. 

One possible argument against microlensing could be 
that, since the two galaxies are fairly distant in projec- 
tion from the quasar pair (~ 35 kpc), there may 
be no stars to account for microlensing. We note that 
evidence of microlensing has previously been fou nd in 
both SDSS J1004+4112 ( e.g.. IRichards et al.l|200l . and 
SDSS J1029+2623 (e.g., IQguri et al.l I2008bj h the two 
large-separation cluster lenses; in those cases however, 
the galaxies believed to host the stars responsible for mi- 
crolensing are located at slightly smaller projected sepa- 
rations of < 20 — 30 h^Q kpc from the quasar images. 

6.3. Intrinsic variability 

16 This does not refer to the increase in A(SDSS)/B and 
A(SDSS)/A at the blue wing of the Mg II BEL, which is due to 
telluric absorption calibrated out in the SDSS spectrum, but not 
in the follow-up spectra. 



We estimate the magnitude of intrin sic quasar vari- 
abilit y using the structure function fe.g.. lYonehara et"al] 
120081 and reference therein) 

/At \ 03 

V = (1 + 0.024M*) — — mag. (3) 
V ^RF / 

This predicts the magnitude change for a quasar of i band 
absolute magnitude Mi, at a given rest frame wavelength 
Arf in units of A, for a rest frame time interval between 
observations of At^p (in days). The structure function 
has been m easured from a quasar sam ple with Mi ~ — 21 
to ~ —30 (jVanden Berk et all 12004) . in the rest frame 
optical/UV, and on a rest frame time range of up to 2 
years. We also note that the characteristic time scale 
for optical variabili ty of quasars is of the order 1 year 
(jlvezic et a l. 2004b), which is comparable with the time 
delay in the source rest frame (0.4 - 2.8 years), in our 
best lensing models. 

We use the structure function to check if it can repro- 
duce spectral changes in quasar A of a factor of ~ 2 in 
flux (~ 0.75 mag), at - 4500 A (~ 1800 A rest frame), 
in the 13 months (~ 160 days rest frame) between the 
original SDSS and the follow-up spectroscopy. Using the 
image magnifications we obtain from the lensing models 
(Table [7]), and the fact that the source rest frame i band 
corresponds to the observed H - Ks bands, we estimate 
the source quasar absolute brightness to be M, ~ —22.0 
to ~ —25.7, depending on the choice of lensing model. 
Since the structure function is larger for fainter quasars, 
we use the former value and obtain an upper estimate 
of Am ~ 0.23 (because this is a magnitude difference, it 
has the same value in the rest frame as in the magnified, 
observed- frame). This value is too small to explain the 
spectral changes at short wavelengths. We note however 
that the structure function might be able to explain the 
flux changes of ~ 0.3 mag (Am = 0.25) in both quasars 
in the z band, in ~ 4 years, and it may also introduce 
flux variations of about 0.26 - 0.23 mag at the observed- 
frame J - Ks bands, during the maximum rest frame 
time interval of the quasar sample, ~ 2 years. We point 
out however that the structure function is a statistical 
average over a sample of quasars, therefore brightness 
changes of larger amplitude are not ruled out in individ- 
ual objects (magnitude differences larger that 0.75 occur 
in about 1% of quasars; iVanden Berk et al.l (f200l ). 

6.4. Estimates of the intrinsic flux ratio 

Our analysis shows that the spectral shapes of A & B 
(in the continua) are compatible with the gravitational 
lensing hypothesis, where the differences are caused by 
differential extinction, microlensing and intrinsic vari- 
ability. Microlensing and extinction are more important 
at short wavelengths than intrinsic variability, and slit 
losses may also contribute to the changes in the spectrum 
of A. Regarding the intrinsic flux ratio, all three spectral 
flux ratios of the two quasars, which we plotted in Figure 
[21 are fairly constant in the range 7200 — 8400 A , with 
A(SDSS)/A (the original SDSS spectrum of A divided by 
the follow-up spectrum of A) fairly close to 1, and A/B, 
A(SDSS)/B ~ 1.4. This latter value is in excellent agree- 
ment with the flux ratio we measure in the J, H and Ks 
bands, where we expect differential extinction and mi- 
crolensing to be weakest. Together, these results suggest 
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that the intrinsic flux ratio, which would ideally be mea- 
sured in the far-infrared or radio, is ~ 1.4 ± 0.2, and all 
three flux-altering physical factors considered are weak 
in the 7200 — 8400 A spectral range. Here we chose the 
error bar on the flux ratio somewhat arbitrarily, slightly 
larger than the scatter in the flux ratio at 7200 — 8400 
A as well as the J, H and K s bands. 

The two quasars were also d etected in the mid-infrared 
by WISE (Wright ct al. 2010), at almost the same epoch 
as the Subaru MOIRCS J, H and Ks observations. The 
flux ratio at 3.4, 4.6 and 12 //m is ~ 1 (Table [g]), but 
we consider this to be broadly consistent with the value 
above, because the quasars were detected at quite low 
S/N ~ 10 — 20, and the fluxes were measured by de- 
blended profile fits, at very poor angular resolution (6'.'1, 
6'.'4, 6'.'5 and 12'/0 in the four bands, respectively). There 
is also the possibility of unknown flux contamination 
from the two bright galaxies. The quasar pair also ap- 
pears to be radio quiet, based on non-detection in the 
Faint Images of the Radio Sky at Twenty-cm survey, 
which is sensitive to > 1 mJy/beam. 

In view of the WISE data, we also consider the pos- 
sibility that the intrinsic flux ratio is ~ 1.0 ± 0.2 (here 
we have broadened the scatter in the values measured at 
at 3.4, 4.6 and 12 /mi). This requires that the flux ratio 
of ~ 1.4 in the continuum, found in the data redward 
of the Mg II BEL, is due to microlensing. This assump- 
tion is still in agreement with the results in ^6.2[ and 
in particular Equation [TJ which is an order-of- magnitude 
estimate. In this interpretation, the flux change at long 
wavelengths occurs on timescales larger than those sam- 
pled by our observations, whereas those at shorter wave- 
lengths are noticeable. This assumption is testable via 
deeper, high resolution mid-infrared observations, as the 
mid-infrared emitting region is expected to be too large 
for microlensing to occur. This interpretation would also 
explain why the flux ratio decreases slightly with increas- 
ing wavelength between the H and Ks bands, reaching 
unity in the WISE bands. 

7. ENVIRONMENT OF THE LENS 

Since our attempts to model the system as a two- image 
gravitationally lensed quasar require a significant addi- 
tional dark matter component, we next look at the en- 
vironment of the lens, in order to search for clues to the 
existence of a possible galaxy cluster. 

We proceed by making a color-color diagram of the 
galaxies in the MOIRCS field of view (~ 1900 x 1660 
kpc at z = 0.899, centered on SDSS J1320+1644), using 
the J, H and Ks bands. We aim to isolate the galaxies 
that are similar in color to Gl and G2, and study their 
spat ial distribution around t he target. We use SExtrac- 
tor (|Bertin fc Arnoutsl [l996') to catalogue and match all 
the objects in the field of view, and we isolate the galax- 
ies by requiring the SExtractor stellarity parameter to 
be less than 0.2 in all three bands. We eliminate a few 
galaxies that are very close together, so that aperture 
photometry would fail, which leads to a sample of 98 
galaxies. 

From Table |H using the aperture photometry, we de- 
termine the colors {J-H,H- Ks) = (0.87 ± 0.07, 0.88 
± 0.05) for Gl, and (0.87 ± 0.08, 0.80 ± 0.05) for G2. 
In color-color space, Gl and G2 lie close to each other, 
and close to the concentration at (J — H, H — Ks) ~ 



(0.9, 0.7). Thus there may be an excess of galaxies at 
the redshift of Gl and G2, indicative of a cluster. 

We attempt a color cut at 0.8 < J — H < 1.1 and 
0.6 < H — Ks < 0.9, which includes the concentration 
of galaxies as well as Gl and G2. In Figure [5] (right), 
we plot the positions of the galaxies selected by this 
color-cut from the MOIRCS field of view centered on 
SDSS J1320+1644; they seem to spatially surround the 
target. We caution however that the broad color cut 
may contain galaxies at different redshift. A multidi- 
mensional color plot would help break the degeneracies, 
but the UH88 bands, even the deep z band, detect only 
a few of the galaxies visible in the MOIRCS bands and 
do not allow a systematic treatment. 

Our analysis therefore hints that Gl and G2 are em- 
bedded at the center of a galaxy cluster, which in turn 
supports the existence of a heavy dark matter halo asso- 
ciated with the cluster, required by our models in fj5] In 
£|Bl we provide an estimate of the content of dark matter. 

8. SUMMARY AND CONCLUSIONS 

We have reported the discovery as part of the SQLS of 
SDSS J1320+1644, a pair of quasars at z = 1.487 with 
two galaxies in between, in a cross-like configuration. We 
have estimated the redshift of the two galaxies as z — 
0.899, based on absorption lines observed in the quasar 
spectra and consistency with the photometric redshifts 
of the galaxies. 

We used the observational constraints to investigate 
whether the system is a two-image gravitationally lensed 
quasar or a binary quasar. We showed that the key to dif- 
ferentiate between the gravitational lens and binary pair 
hypotheses is whether or not there exists a dark matter 
halo at the location of the system, capable of boosting the 
gravitational potential of the galaxies. We estimated the 
required virial mass or velocity dispersion of this halo, 
as well as the expected time delays and image magni- 
fications, using a variety of SIS and NFW profiles. We 
were also able to show that the halo can only be centered 
on the brighter, elliptical galaxy Gl or somewhere in be- 
tween the two quasars (but not on the disk-like galaxy 
G2), as well as to set constraints on the velocity dis- 
persions of these two galaxies, and an additional, fainter 
third galaxy. The third galaxy may be located at ei- 
ther z = 0.899 or z = 1.168. In the first case, it would 
further support the existence of a halo indicative of a 
galaxy group or cluster, and in the second it would ex- 
plain the presence of a strong absorption system in one 
of the quasars. 

We have found a substantial number of galaxies in the 
field, surrounding the quasar pair, with colors consistent 
with z = 0.899 (based however on only three filters). 
If this is a galaxy cluster, the quasar pair would reside 
behind the center of the cluster potential, which would 
provide the required dark matter halo lens. The calcu- 
lated mass-to-light ratio assuming the lensing hypothesis 
( £|Bl) is in agreement with what would be expected from 
a cluster-scale lens. 

We have performed a comparative analysis of the spec- 
tra of the two quasars (as well as inferred the intrinsic 
flux ratio of the two quasars), and concluded that all 
observed differences can be attributed to a combination 
of extinction, microlensing and intrinsic variability, in 
support of the gravitational lensing hypothesis, with the 



microlensing being dominant. 

The most similar known system to our target is ar- 
guably Q0957+561 (jWalsh et all 119791) . which has two 
quasar images separated by 6'.' 17 and one lensing galaxy 
that is part of a confirmed cluster that boosts the image 
separation. That system shows, in addition, evidence of 
strongly lensed arcs, which ar e the extended image s of 
the lensed quasar host galaxy ()Bernstein et al.lll993D . 

We conclude that SDSS J1320+1644 is a probable 
gravitationally lensed quasar, although we are unable to 
prove this beyond doubt, based on our available data. 
We propose the following observations as a means of es- 
tablishing the true nature of this object unambiguously: 

1. As the most time-effective way of testing the lens- 
ing hypothesis, we suggest deep adaptive optics 
or Hubble Space Telescope observations to study 
the morphology of the quasar host galaxies, and 
look for extended arc-like features, characteristic 
of gravitational lensing. We note that there is a 
16.6 mag (R band) star ~ 58'.' 4 from the target, 
which could be used as a tip-tilt star for the Laser 
Guide Star Adaptive Optics capabilities of Subaru 
Telescope. 

2. Deep, wide field optical imaging of the region could 
be used to detect weak lensing associated with the 
possible lensing cluster. 

3. Deep X-ray observations would test the presence 
of the galaxy group/cluster by providing an inde- 
pendent measure of the velocity dispersion of intr- 
acluster g well as determining its dynamical 
center. 

4. Spectroscopy of the two brightest galaxies would 
determine if they are indeed at the same redshift, 
consistent with the spectral absorption lines, as 
assumed in our galaxy cluster hypothesis. These 
galaxies are faint in the V, R, I, z and SDSS bands, 
and therefore the photometric redshifts are fairly 
uncertain. Multi-object spectroscopy would pro- 
vide a means to check the redshift of the surround- 
ing galaxies, and therefore either consistently prove 
or disprove the cluster hypothesis. At the very 
least, deeper imaging in other bands is required to 
constrain the photometric redshift of the surround- 
ing galaxies. 

5. Long-time monitoring of the quasar fluxes to look 
for correlated variations (a time delay) would con- 
vincingly reveal the nature of this system, as long 
as the correlated variations are not suppressed by 
microlensing. Although we do not possess suffi- 
cient constraints to estimate a robust time delay, 
our best estimates are ^1 — 7 years, based on our 
lensing models. 

Should this system indeed prove to be a gravitational 
lens, it would be the third-largest separation gravitation- 
ally lensed quasar in the SQLS, and the largest separa- 
tion two-image lensed quasar known. These large sep- 
aration lenses are important to include in a complete 
statistical sample s uch as the one provided by the SQLS 
(|Inada et alj 120121) . in order to constrain the hierarchi- 
cal structure formation at cluster mass scales. They are 
also exceedingly rare, because the probability of quasars 
strongly lensed by clusters is 1-2 orders of magnitudes 



smaller than tha t of quasars lensed by galaxies (e.g., 
Ilnada et all [20121 ) . 
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Table 1 

Original SDSS photometry 



Object & designation 


u 


9 


r 


i 


z 


SDSS J132059. 17+164402.59 (A) 


20.16+0.04 
20.10+0.06 


19.58+0.02 
19.54+0.02 


19.16+0.02 
19.17+0.02 


18.88+0.02 
18.84+0.02 


18.86+0.04 
18.84+0.06 


SDSS J 132059.73+164405.6 (B) 


18.93+0.02 
18.91+0.02 


18.84+0.02 
18.78+0.01 


18.74+0.02 
18.74+0.02 


18.62+0.02 
18.57+0.02 


18.67+0.03 
18.69+0.06 


SDSS J132059.37+164406.66 (G2) 


25.05+0.79 
24.53+2.84 


24.13+0.36 
23.99+0.58 


23.01+0.22 
23.08+0.41 


22.10+0.17 
21.96+0.27 


21.77+0.47 
21.60+0.80 



Note. — The AB magnitudes are c orrec ted for Galactic extinction following Schlcgcl ct al. (1998). 
The object designations are those from i|3,2l The first line for each object gives the PSF magnitudes (A, 
B), or best model magnitudes (G2), and the second line gives aperture photometry magnitudes. Objects 
G3 and G5 (see Figure [4j were also identified as photometric targets by the SDSS pipeline, but their 
magnitudes are not reproduced, as they are considered unreliable. 



Table 2 

Summary of follow-up imaging observations 



Filter 


Instrument 


Exposure 


Airmass 


Seeing 


Zero-point uncertainty 


Observation date 


V 


UH88/Tek2k 


3 X 100 s 


1.20 


1.0 - 1.2 


0.10 


2009 April 15 


R 


UH88/Tek2k 


4 x 100 s 


1.40 


1.0 - 1.2 


0.08 


2009 April 15 


I 


UH88/Tek2k 


7 x 100 s 


1.20 


1.0 - 1.2 


0.07 


2009 April 15 


z 


UH88/Tek2k 


8 x 100 s 


1.20 


1.0 - 1.2 


0.07 


2009 April 15 


z 


UH88/WFGS2 


24 x 180 s 


1.14- 1.53 


1.0 - 1.2 


0.04 


2009 April 17 


.) 


Subaru/MOIRCS 


4 x 150 s 


1.03 


0.4 


0.01 


2010 April 2 


H 


Subaru/MOIRCS 


4 x 120 s 


1.05 


0.1 


0.01 


2010 April 2 


Ks 


Subaru/MOIRCS 


4 x 150 s 


1.09 


0.4 


0.01 


2010 April 2 


Note. 


— Zero-point uncertainties in the V, R, I and 


z filters include the scatter between the values obtained 



from the 10 stars used in the estimation, as well as the Lupton (2005) transformations. In the J, H and Ks 
filters, the uncertainties include those quoted in the catalogue and the aperture photometry errors. There the 
actual error may be ~ 0.05 mag, if the differences between catalogues and the values obtained from different 
frames of the standard star are included. Among the two z band observations, only the deeper one is considered 
throughout this article. 



Table 3 

Astrometry results for SDSS J1320+1644 



Object 


Ax [arcsec] 


Ay [arcsec] 


A 


0.000 ± 0.002 


0.000 ± 0.002 


B 


-8.066 ± 0.002 


2.939 ± 0.002 


Gl 


-4.991 ± 0.006 


0.117 ± 0.006 


G2 


-2.960 ± 0.006 


3.843 ± 0.007 


Gl 


-9.169 ± 0.019 


5.173 ± 0.024 


G3 


-1.682 ± 0.018 


-5.166 ± 0.022 


G5 


-12.177 ± 0.021 


7.176 ± 0.020 


Note. 


— Astrometry 


of the SDSS 



J1320+1644 system, determined in the 
MOIRCS Ks band (the Sersic index of the 
galaxies has been left unconstrained in the 
GALFIT modeling). The positive directions 
of x and y are defined towards West and 
North, respectively. The quoted errors are the 
GALFIT statistical errors (A, B, Gl, G2), and 
those determined with the IRAF PHOT task 
(G3. G4, G5). 
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Figure 1. SDSS gri color-composite (Lupton ct al. 2004) centered on the SDSS J1320+1644 pair, obtained on 2005, June 6. The image 
scale is 50" X 50". North is up and East is to the left. The exposure was 54 s in each band and the SDSS pixel scale is ty/395. The two 
quasars are the blue objects near the center. [See the electronic edition of the Journal for a color version of this figure.] 




4000 5000 6000 7000 8000 9000 

^obs [A] 

Figure 2. Top: Spectra of the two stellar components (quasars) in SDSS J1320+1644, A and B. For component A, we plot the original 
SDSS spectrum (blue), as well as the ARC 3.5 m follow-up spectrum (black). The quasar designation corresponds to Figure [4] (object A is 
the component to the South- West, in Figure[TJ. The emission and absorption lines we identified are marked. Two absorption systems in 
the follow-up spectra, at the blue wing of Mg II and at ~ 7600 A, are due to atmospheric absorption (these have been calibrated out of the 
original SDSS spectrum). Bottom: the ratio of all three spectral pairs. 
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Figure 4. Subaru Telescope imaging of SDSS J1320+1644. The image scale (without the close-up) and orientation correspond to Figure 
[T] This is a small fraction of the MOIRCS field of view centered on SD SS J1320+1644 (see the panel to the right in Figure [Sj . The image 
is a color-composite of the J, H and Ks frames, using the algorithm of Lupton ct al. (2004). [See the electronic edition of the Journal for 
a color version of this figure.] 
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Table 4 

Follow-up photometry of the SDSS J1320+1644 system 



Object 


V 


R 


1 


z 


J 


H 


K s 


A 


19.23 ± 0.01 
19.20 ± 0.01 


18.68 ± 0.01 
18.71 ± 0.01 


18.33 ± 0.01 
18.31 ± 0.01 


18.62 ± 0.01 
18.71 ± 0.01 


17.58 ± 0.01 
17.62 ± 0.02 


16.79 ± 0.01 
16.86 ± 0.01 


16.62 ± 0.01 
16.67 ± 0.01 


B 


19.30 ± 0.01 
19.28 ± 0.01 


18.92 ± 0.01 
18.94 ± 0.01 


18.67 ± 0.01 
18.66 ± 0.01 


19.00 ± 0.01 
19.10 ± 0.02 


18.05 ± 0.01 
18.04 ± 0.02 


17.29 ± 0.01 
17.29 ± 0.02 


17.03 ± 0.01 
16.98 ± 0.02 


Gl 


23.63 ± 0.37 
23.02 ± 0.23 


23.01 ± 0.22 
22.85 ± 0.16 


21.66 ± 0.11 
21.47 ± 0.09 


21.62 ± 0.10 
21.83 ± 0.12 


19.96 ± 0.02 
20.00 ± 0.06 


18.83 ± 0.11 
19.13 ± 0.04 


18.13 ± 0.04 
18.25 ± 0.03 


G2 


23.16 ± 0.24 
23.25 ± 0.28 


23.04 ± 0.22 
23.52 ± 0.29 


22.12 ± 0.16 
22.52 ± 0.23 


21.51 ± 0.08 
21.65 ± 0.10 


20.14 ± 0.02 
20.20 ± 0.06 


19.24 ± 0.02 
19.33 ± 0.05 


18.50 ± 0.01 
18.54 ± 0.04 


G3 










21.42 ± 0.12 


20.24 ± 0.10 


20.07 ± 0.10 


G4 










21.86 ± 0.14 


21.05 ± 0.14 


20.35 ± 0.10 


G5 










22.23 ± 0.17 


21.08 ± 0.14 


20.80 ± 0.13 



Note. — The values are in magnitudes (all in the Vega system, except for the z band, in the AB system). For 
A, B, Gl and G2, the first line indicates model magnitudes provided by GALFIT (morphological fit), and the second 
line aperture photometry obtain ed with the IRAF P HOT task. For G3, G4 and G5, only aperture photometry results 
are given. Galactic extinction (Schlcgcl ct al. 1998) and atmospheric extinction have been corrected. Quoted errors 
are statistical errors only. Errors do not include uncertainties in the photometric zero-points (Table [2jl, in the PSF 
or due to the galaxy modeling parameters. Aperture photometry has been measured in all bands with a fixed ~ 3'.'1 
diameter aperture (~ 1"9 for G4 and G5), selected to maximize the amount of light from the target, as well as 
minimize contamination from the other objects. In the V , R, I and z bands, contamination may be higher due to the 
larger seeing. In order for the models to converge, the morphological parameters of the Sersic profiles fitted to these 
four bands were constrained to the values determined in the Ks band (i.e., Sersic index fixed to the closest canonical 
value of 1 or 4, effective radius). 



Table 5 

Morphological parameters of the two bright galaxies 



Object 


Effective radius [arcsec] 


Sersic index 


Axis ratio (b/a) 


Position angle [deg] 


Gl 


0.24±0.01 


4.72±0.76 


0.60±0.04 


-3.6±3.5 


G2 


0.37±0.01 


0.57±0.07 


0.66±0.02 


2.6±2.7 



Note. — Morphological parameters are fitted by GALFIT in the Ks band, with attached 
statistical errors. The position angle is measured from North towards East. 




Figure 5. Photometric redshift probability distributions, normalized to unit area, for the two galaxies Gl (left) and G2 (right). The 
probability distributions are calculated with EAzY using the magnitudes in Tables[l]and|4](see also J4]for details) for the elliptical template, 
Sbc template, or marginalized over all spectral templates. The redshift of the quasar pair, as well as those of the absorption line systems 
identified in the quasar spectra, are marked with vertical lines. 
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Table 6 

Best-fit photometric redshifts 







JJCOv 111 "6 


lcr 


. . 

in 1 its 


2cr limits 


X /d.o.f. l a -°- r -; 


Gl, E (EAzY) 


VRIzJHKs 


0.81 


0.67 


- 0.96 


0.53 - 


1.10 


0.34 (6) 


Gl, E (HyperZ) 


VRIzJHKs 


0.87 


0.77 


- 0.94 


0.71 - 


0.98 


0.65 (6) 


Gl, E (EAzY) 


JHKs 


0.83 


0.65 


- 1.02 


0.38 - 


1.20 


0.17 (2) 


Gl, E (HyperZ) 


JHKs 


0.93 


0.79 


- 1.04 


0.71 - 


1.19 


0.03 (2) 


G2, Sbc (EAzY) 


Mgri VRIzJHKs 


0.80 


0.71 


- 0.89 


0.64 - 


0.97 


0.46 (10) 


G2, Sbc (HyperZ) 


ugriVRIzJHKs 


0.86 


0.80 


- 0.91 


0.67 - 


0.95 


1.97 (10) 


G2, Sbc (EAzY) 


JHKs 


0.79 


0.57 


- 0.97 


0.38 - 


1.32 


0.29 (2) 


G4, all (EAzY) 


JHKs 


1.11 


0.56 


- 1.58 


0.29 - 


1.98 


0.05 (2) 


G3, all (EAzY) 


JHKs 


0.55 


0.47 


- 1.10 


0.27 - 


1.55 


2.07 (2) 


G5, all (EAzY) 


JHKs 


0.60 


0.52 


- 1.45 


0.30 - 


1.94 


0.84 (2) 




Figure 6. Gravitational lensing models. Top left: A and B are two images of a source lensed primarily by a large velocity dispersion 
SIS (with shear) profile centered on Gl; top right: same as previous, but the position of the primary SIS (without shear) lens is not fixed; 
bottom left: A and B are two images of a source lensed primarily by a large elliptical NFW profile centered on Gl; bottom right: same as 
previous, but the position of the spherical NFW lens is not fixed. Critical lines and caustics are drawn in blue (outer curves) and light 
green (inner curves), respectively, with the velocity dispersions of Gl, G2 and G4 fixed at their most probable values (237, 163 and 118 km 
s —1 , respectively; see [[5] and unless the adopted model requires otherwise. Only NFW profiles with c„; r = 6 are drawn. Symbols 

denote image, source and lens positions (red squares, magenta triangles and black crosses, respectively). The black regions close to the 
center of the large caustic of the free position models show the locations of this primary lens, as the velocity dispersions of Gl, G2 and G4 
are varied over the allowed lcr range (x 2 < 2.3). [See the electronic edition of the Journal for a color version of this figure.] 
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Table 7 

Parameters of the best-fit lensing models 



model 


a [km s _1 ] or M [/i 7o 'M ] 




e or 7 


8 e or 0-y [deg] 




At [days] 


SIS+7 at Gl 


710 ±8 km s" 1 




0.026 ±0.011 


20 ±20 


18 ±2 


911 ±21 


SIS free 


645 ±25 km s" 1 








37 ±29 


-860 ±460 / 
1630 ± 940 


elliptical NFW at Gl 


(5.5 ± 0.5) x 10 13 h^Q Mq (c vir 
(4.2 ±0.7) x 10 14 h^M @ (c vir 
(3.15 ±0.35) x 10 15 hjnM® (c vl 


= 30) 
= 6) 
r = 2) 


0.03 ±0.01 
0.018 ±0.012 
0.01 ±0.01 


20 ±13 
~ - 180 
~ - 180 


13.5 ±0.5 
60 ± 12 
147 ± 43 


1095 ± 15 
510 ±40 
320 ± 41 


NFW free 


(4.9 ±0.4) x 10 13 /ifo M© ( c vir 
(4.0 ±0.4) x 10 14 h?}M @ (c vir 
(3.2 ±0.1) x 10 15 hfo Af (c vir 


= 30) 
= 6) 
= 2) 






41 ± 18 
99 ±46 
175 ± 75 


-638 ± 262 
-547 ± 223 
-380 ± 80 



Note. — The range on each parameter is obtained by changing the velocity dispersions of Gl, G2 and G4 over the allowed 
la range (see ij5j. e and 7 are the ellipticity and shear, respectively, and 8 e , 8-y arc the respective position angles (measured 
East of North), /itot is the total magnification of images A and B, and the predicted time delay At is positive if image A leads. 



cq 



A 2005/6/6 SDSS 
B 2005/6/6 SDSS 
A 2009/4/15,17 UH88 
B 2009/4/15,17 UH88 
A 2010/4/2 Subaru 
B 2010/4/2 Subaru 




Figure 7. Broadband photometric flux and flux ratios of quasars A and B, determined from the aperture photometry in Tables [T] and 
[4] which include curve of growth corrections. The photometric bands are marked, and the horizontal bars of the flux ratios indicate the 
width of the filters. We use the photometric errors corr espo nding to the zero-point uncertainties, from Table [2] The WISE mid-infrared 
data from Table [8] which is considered less reliable (see £|6,3l L is not plotted. [See the electronic edition of the Journal for a color version 
of this figure.] 



Table 8 

WISE mid-infrared photometry of the SDSS J1320±1644 quasar pair 



Object & flux ratio 


Wl (3.4 fim) 


W2 (4.6 ftm) 


W3 (12 urn) 


W4 (22 /im) 


A 


15.32 ± 0.05 


13.99 ±0.05 


10.92 ±0.10 


9.06 ±0.53 


B 


15.32 ±0.05 


14.04 ± 0.05 


10.85 ±0.09 


7.97 ±0.20 


A/B 


1.00 ±0.07 


1.05 ±0.07 


0.94 ±0.12 


0.53 ±0.19 



Note. — The deblended profile-fit magnitudes are in the Vega system. The 
observations were performed between 2010, January 7 and 2010, June 30. The 
values in W4 are particularly unreliable, as the quasars are detected at S/N ~ 2 — 5, 
with an angular resolution of 12'.'0, larger than the quasar separation. 
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Figure 8. Color-color (left) and spatial distribution diagrams (right) of the galaxies in the Subaru MOIRCS field of view. The colors and 
associated error bars are based on aperture photometry (MAG_APERT) determined with SExtractor. The color-cut region used to isolate 
galaxies possibly associated with the candidate cluster at z ~ 0.9 is marked with dashed lines. The spatial distribution diagram identifies 
the galaxies whose colors match the color-cut region. We use the field of view of only one of the MOIRCS detectors, centered on the target. 
East is up and North is to the right. 100 arcsec at z = 0.9 corresponds to a proper scale of 790 h^ kpc. [See the electronic edition of the 
Journal for a color version of this figure.] 
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APPENDIX 

A. ADDITIONAL CONSIDERATIONS ON THE LENS MODELING 
A.l. A lens model for G4 

By considering that G4 is located at z = 1.168 and applying the F-J and T-F laws, the velocity dispersion of 
this galaxy is a = 175 ± 40 km s" 1 or 125 ± 12 km s _1 , respectively. Assuming that the effect of G4 on the 
lensing configuration is small, we can treat G4 as an SIS lens situated at z = 0.899, but with a smaller effective 
velocity dispersion. This eliminates the need to introduce multiple lens planes and is justified by the fact that the 
convergence from G4 at B (arbitrarily taken as the object closest to G4, but projected to z — 1.168) is as small as 
k ~ 6 Ein /2d B -G4 ~ 0.028 ± 0.012 (F-J) or - 0.014 ± 0.002 (T-F), where 6> E i n is the Einstein radius of the lens G4 
estimated from the velocit y dispersion (assumin g an SIS profile) and cLb-ga — 2'.'49 is the projected angular distance 
between B and G4. From [Keeton e t al. (2003), using n (convergence) = 7 (shear) for the SIS profile, the effective 
convergence of G4 at z = 0.899 is K eB = (1 - P)k/(1 - 2fin) ~ 0.012 ± 0.005 (F-J), or - 0.006 ± 0.002 (T-F), where 
fi = -Di2-Dos/(-Do2-Dis), Dij = D{z il Zj) are angular diameter distances, z\ = 0.899, z 2 = 1.168, and O, S refer to the 
observer and source, respectively. The effective convergence corresponds to an SIS profile at G4 with effective velocity 
dispersion ~ 81 ± 19 km s _1 (F-J) or ~ 57 ± 5 km s _1 (T-F). We conclude that, considering all possibilities, an SIS 
model at z = 0.899 for G4 would have a velocity dispersion in the la range of ~ 118 ± 66 km s _1 . 



A. 2. Dependence on c v i r 

Here we investigate how the choice of the concentration parameter c V i r of the NFW profiles affects our results in 
We explore c V i r in the range from 2 to 30, and we optimize all other model parameters. Increasing c V i r has the effect of 
decreasing the total magnification of the images, and increasing the time delay (Table [7]). We plot the resulting c vir - 
virial mass relation in Figure [9l We obtain virtually the same degeneracy for the free NFW profile as for the one fixed 
at G2. For comparison, we also plot some recently published results on the c V i r - virial mass relation: t he results of 
N— bod y simulations bv lDuffv et al.l (|2008l ). as well as the theoretical and observational results in ferred bvlQguri et al.l 
(|2012aT> from fitting lensing clusters at z ~ 0.45. The strong lens selected sample behind the lOguri et al.l ( 2012al) 



relations is known to be biased towards overestimating the concentrations due to halo triaxiality. If SDSS J1320+1644 
is indeed a lens system, it should include the same bias, as it is selected as a strong lens. All these relations are 
observed to intersect the SDSS J1320+1644 curves. 



A. 3. Caveats on the velocity dispersions of the galaxies 

We discuss two caveats. A first caveat is that, when estimating the range of velocity dispersion for each galaxy in <j5j 
we have not taken into account the p ossible evolut i on of the T-F and F-J with redshift . Althoug h there is no c onsen sus 
in the literature in this regard (e.g., iBohm et all (|2004h & iFernandez Lorenzo! (|2010f ) for T-F, iRusin et all (|2003| ) & 
iTreu et all (|2006l ) for F-J), the tendency is that galaxies had smaller rotational velocities and velocity dispersions in 
the past for the same luminosity. This would reduce even more the contribution of the galaxies to the mass models 
dominated by the dark matter halo. 

A second caveat applies to G4, in the case that this galaxy is located at z = 1.168. In this case, G4 would also 
be lensed by the dark matter halo at z = 0.899, and therefore magnified by a factor ~ 1.5 — 3, depending on the 
assumed lensing model. This would mean that the galaxy is intrinsically fainter, and its contribution to the mass 
models would be further reduced. On the other hand, the undeflected position of G4 would be closer to the main 
deflector, increasing its influence. We have checked that the best-fit models considered in this section (from which 
we exclude the contribution of G4) would not produce multiple images of this galaxy (z — 1.168 being closer to the 
redshift of the main lens than z = 1.487, the critical surface density for multiple images is increased). 

B. MASS-TO-LIGHT RATIOS 

Here we compute the mass-to-light (M/L) ratio by assuming that SDSS J1320+1644 is a lensed system, and compare 
it with the M/L ratio of the other lenses discovered by the SQLSFI ()Inada et al.ll2012l ). We fit each lens with an SIS 
profile, in which case the Einstein radius 9e is simply half the image separation. For the quasars with more than two 
images, we consider 9e to be half the largest separation between the images. The mass inside the Einstein radius is 
simply determined as 

M « e E ) = B°^L?3l „ 1.23 x 10 8 h^u @ ( D °^ s/DhS ) (V^-V. (Bl) 

V ; D LS AG 70 V 1 /ifo 1 Mpc J Vlarcsecy v ' 

To determine the luminosity of the SQLS lenses, we referred to the discovery papers (jlnada et al.ll2012l and references 
therein) and estimated the rest frame R band magnitudes of all the galaxies with available photometry, located 
between the images, using the HyperZ algorithm. To convert the rest frame magnitudes into luminosities we used 

17 http : / /www-utap .phys . s . u-tokyo . ac . jp/~ sdss/sqls/lens . html 
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Mr = 4.48 — 2.5 log L/Lq^r, where M Qj r = 4.48 is the R band magnitude of the Sun. The M/L ratios were computed 
in terms of the solar M/L ratio, T© = 1 M Q /L Q . 

We attached no error bars to the results of these rough calculations. Not accounting for the effects of lens ellip- 
ticity and external shear for each (mostly double- image) system should introduce an error of ~ 10% on the masses 
(R usin fc K ochanck 2005]). Also, the redshifts of several lenses are not accurately known, which influences both the 
masses and the inferred luminosities. Photometry for most lenses is only available in three bands (we eliminate 10 
objects for which the lens is detected in only one band), and their morphology is typically unknown. It is therefore 
difficult to choose the appropriate CWW spectral template for calculating the rest frame R band l uminosity. Howeve r 
we know that the majority of lens galaxies are ellipticals, as these have larger velocity dispersions (jTurner et al.lll98~4T) . 
and so wc generally assume the elliptical template, unless the lens is known to be a spiral, or the fit of another tem- 
plate is much better. Also, luminosities should be more accurate for lenses at z ~ 0.2 - 0.3, because at this redshift 
the rest frame R band is rcdshiftcd into the observed-frame I band, where observations are usually available. For 
SDSS J1320+1644, the lens redshift estimate means that the rest frame R band is redshifted into the J band, where 
photometry is available. 

The results of our M/L estimates are shown in Figure [101 in which we plot the mass-to- light ratios against the 
reduced image separations. The reduced image separation A# rc d = A8 ■ Dos/D^g scales out the dependence on the 
source and lens distance, thus representing a physical property of the lensing object. 

The image separation is the most important observable, because it reflects the depth of the gravitational potential 
of the lensing object and therefore the structure responsible for the lensing phenomenon. Since clusters are dominated 
by dark matter much more than normal galaxies, we expect that, as the image separation increases from the galaxy 
lens to the cluster lens range, so does the M/L ratio. Accordingly, the result we obtain in Figure [T0l shows a positive 
correlation between the image separation and M/L. Most small-separation lenses at A9 ~ 1" are clustered at M/L 
~ 1 — 10 T(p, consistent w ith the estimates from the dynamics of elliptical galaxies and their stellar populations 
(jPadmanabhan et al.ll200l . whereas the two known cluster-scale lenses SDSS J1004+4112 and SDSS J1029+2623 
both have a larger M/L. We obtain a similarly large M/L for SDSS J1320+1644. 



A quasar - galaxy cross 
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Figure 9. The c v i r - virial mass degeneracy of the NFW profiles. The thick black lines show publishe d re l ations from both simulations 
and observations, while the thin lines show the corresponding la ranges. A caveat is that the Oguri ct al. (2012a) theoretical curve has 
been derived for a cluster redshift z = 0.45, which corresponds to the mean redshift of their observed cluster sample. The blue and red 
curves show the results of these work, for the NFW profile fixed at Gl, as well as with unconstrained (free) location, respectively. [See the 
electronic edition of the Journal for a color version of this figure.] 
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Figure 10. Mass-to-light ratios in the rest frame R band for 52 SQLS gravitationally lensed quasars, as well as for SDSS J1320+1644 
(assuming it is also a lensed quasar). Here A0 re£ ; is the reduced image separation (see JB}. 



